Impact of supernova dynamics on the vp-process 
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ABSTRACT 

We study the impact of the late time dynamical evolution of ejecta from core-collapse supernovae 
on y/7-process nucleosynthesis. Our results are based on hydrodynamical simulations of neutrino wind 
ejecta. Motivated by recent two-dimensional wind simulations, we vary the dynamical evolution during 
the y/7-process and show that final abundances strongly depend on the temperature evolution. When the 
expansion is very fast, there is not enough time for antineutrino absorption on protons to produce enough 
neutrons to overcome the -decay waiting points and no heavy elements beyond A = 64 are produced. 
The wind termination shock or reverse shock dramatically reduces the expansion speed of the ejecta. This 
extends the period during which matter remains at relatively high temperatures and is exposed to high 
neutrino fluxes, thus allowing for further (/?, y) and (n, p) reactions to occur and to synthesize elements 
beyond iron. We find that the y;?-process starts to efficiently produce heavy elements only when the 
temperature drops below ~ 3 GK. At higher temperatures, due to the low alpha separation energy of ^^Zn 
{Sa = 2.7 MeV) the reaction ^^Cu(;?, CK)^^Ni is faster than the reaction ^^Cu(j9, y)^^Zx\. This results in the 
closed NiCu cycle that we identify and discuss here for the first time. We also investigate the late phase of 
the y/?-process when the temperatures become too low to maintain proton captures. Depending on the late 
neutron density, the evolution to stability is dominated by decays or by (n, y) reactions. In the latter 
case, the matter flow can even reach the neutron-rich side of stability and the isotopic composition of a 
given element is then dominated by neutron-rich isotopes. 

Subject headings: Nuclear reactions, nucleosynthesis, abundances — supernovae: general 

1. Introduction interact with the stellar matter and deposit energy in 

the outer layers of the proto-neutron star leading to 
Neutrino-driven winds from core-collapse super- ^ supersonic outflow known as neutrino-driven wind 

nova explosions contribute to the synthesis of ele- ([Puncan et al. 1986 ). Although neutrino-driven winds 

ments beyond iron. After the explosion, the hot proto- ^gj-e considered the site where heavy element s are 

neutron star cools emitting neutrinos. These neutrinos produced by the r-process ( Wooslev et al. "1994), re- 
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cent simulations dArcones et al .. .2007.:„Hudepohl et"aD 
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2Q1QI: iFischer et al.ll201Ql: [Roberts etal]l201Qh cannot 
reproduce the extreme conditions required for produc- 
ing heavy r-process elements (see e.g., Hoffman et al.l 
[1997!: lOtsuki et al.ll20QQl: [Thompson et al.ll2QQll) . The 
wind entropy is too low (less than lOO^e/nuc) and, 
even more significant, the ejecta is proton rich (the 
elec tron fraction 7^ r e mains ab ove 0.5 during s econd s, 
see iHiidepohl et aP feoioh : iFischer et all tOld) ). 



Even if the r-process does not take place in every 
neutrino-driven wind, lighter heavy elements (e.g., 
Sr, Y, Zr) ca n be syn thesized in t his en vironment as 
suggested by lOian &^W asserbur3 j2001 ). In proton- 
rich conditions iFrohlich et al.l 12006) showe d that ele- 
ments b eyond ^"^Ge can be synthesized. Wanajo et al 
(1201 Ibh found Sr, Y, Zr in small pockets of neutron 



rich material ejected after the explosion of low mass 
progenitors. Recently, lArcones & MontesI (|201l|) 
performed a systematic nucleosynthesis study that 
strongly supports the production of lighter heavy el- 
ements in proton- and neutron-rich neutrino-driven 
winds. 

In proton-rich winds, charged particle reactions (al- 
pha and proton captures) build nuclei up to ^^Ni and 
even up to ^"^Ge once the temperature drops below 
3 GK. Due to their long beta-decay lifetimes and 
low thresholds for proton capture, the nuclei ^^Ni 
and ^"^Ge act as bottlenecks that inhibit the produc- 
tion of heavier elements. In the y/7-process, their de- 
cay is sped up by (n, p) reactions, with the neutrons 
produced by antineutrino absorption on the abundant 
free protons. This allows for the production of el- 
ements beyond iron and may explain the origin of 



light p-nuclei (Frohlich et a l. 2006; Pruet et al. 
Wanajol I2OO6I: IWanajo ^\\ boilah . The synthesis 



of elements by the y/?-process depends thus on neu- 
trino spectra and luminosities but also on the dynam- 
ical evolution as matter expands through the slow, 
early supernova ejecta. This produces a wind ter- 
mination shock or reverse shock where kinetic en- 



ergy is transformed into internal energy jArcones et al 



ergy 
I2OO7 



2007 ). The reverse shock has a big impact on the 
nucleosynthesis because temperature and density in- 
crease and the expansion is strongly decelerated. This 
hydrodynamical feature has been extensively studied 
for r-process nucleosynthesis (Qian & Woosley 1996; 
Sumivoshi et al.ll2000l: IWanaioll20Q7l: IPanov & Jankal 



2009': 'Arcones & Martme z-PinedolEoTlh . Recently, 
Wanajo et al. (201 la) have also explored the relevance 
of the reverse shock on the y/7-process. Motivated by 
their work and by new 2D hydrodynamical simulations 



of the neutrino-driven wind ( Arcones & Jankal I20 1 1 ). 
we investigate here the wind termination shock to gain 
further insights on the dynamical evolution relevant for 
the y/7-process. 

Here we use a trajectory from hydrodynamical sim- 
ulations (Sect. [2]) combined with a complete nucle- 
osynthesis network (Sect.[3l). In Sect.|4]we present our 
results where we analyze the impact of the wind termi- 
nation temperature (Sect. 14.11) . the temperature jump 
at the reverse shock (Sect. 14.21) . and the influence of 
late temperature evolution and the decay to stability 
(^Sect l4.3l) . Our conclusions are summarized in Sect. [51 

2. Long-time dynamical evolution 



Our study is based on one tr ajectory obtain e d from 
hydrodynamical simulations by lArcones et al.l (|2007|) . 
These simulations follow the supernova explosion and 
the subsequent neutrino-driven wind. 
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Fig. 1 . — Evolution of density and temperature along a 
tr ajectory ejected at 2s after bounce of model MlSUrl 
in 'Arco nes et al.l (l2007k . The reverse shock is at T^wt ~ 
0.6 GK and produces a jump in density and tempera- 
ture and a sudden decrease of the expansion velocity. 

Figure [T] presents the evolution of density and tem- 
perature for the chosen trajectory. In the following, 
we always use the same initial evolution and mod- 
ify it at temperatures below 3 GK. We assume the 
wind terminates at a temperature T^t and the evo- 
lution therea fter is varied using the pre scription in- 
troduced by Arcones & Martmez-Pinedo (2011) and 
motivated by ID and 2D hydrodynamical simu lations 



(lArcones et al.l 120071: lArcones & Jankal l201lk . The 



thermodynamical conditions at the termination of the 
wind depend on both the intrinsic properties of the 
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wind, like its velocity and mass outflow rate, and on 
the pressure of the slow moving ejecta. The wind 
properties are determined by the neutrino energies 
and luminositie s and by the neutron star mass and 
radius (see e.g.JOian & Woosleylll996h . These condi- 
tions depend mainly on the nuclear eq uation of state 



and a r e similar for different progenitors (|Arcones et al 
2QQ7I: iFischer et al] I201Q|) . The properties of the 



slow supernova ejecta are related to the progenitor 
and become very anis otropic during the evolution 
dArcones & Jankall201ll) . 

In our calculations, we use a simple model that re- 
produces the main features seen in hydrodynamical 
simulations. Once the wind reaches the early super- 
nova ejecta we use the Rankine-Hugoniot conditions 
to determine the behavior of temperature, density, and 
velocity. When the wind moves supersonically these 
quantities jump at the reverse shock. The evolution 
after such discontinuity is determined assuming con- 
stant density and temperature during a time Ar. As 
the mass outflow is constant (M = 47ir^vp) the ve- 
locity drops as V oc r"^. At later times, the velocity 
stay s constant and the den sity decreases as p oc r"^ 
(see lArcones & Martmez-P inedo (2011) for more de- 
tails). To account for the differences found in the late 
evolution in two-dimensional simulations, we will use 
different values of and explore the impact on the 
nucleosynthesis. 

3. Nucleosynthesis network 

The evolution of the cor nposition is calculate d us- 
ing a full reaction network dprohlich et al .1120061) . All 
calculations start when the temperature drops below 
10 GK and are followed until the temperature reaches 
0.01 GK. We assume the initial composition to be 
determined from nuclear statistical equilibrium for 
a fixed electron fraction Ye = 0.52. The reaction 
network includes 1869 nuclei from free nucleons to 
dysprosium (Z = 66). Neutral and charged parti- 
cle reactions are taken from the REACLIB compi- 
lat ion, containing theoretical stat istical-model rates 
by iRauscher & Thiele mann (2000) (NON-SMO KER) 
and experimental rates by Angulo et al.l (Il999|) . The 
theoretical weak interactions rates are taken from 
Langanke & Martinez-Pinedol (1200 1 *) and F uller et al. 
(Il982h . Where available, experimental beta-decay 



rates are used from NuDat2Q supplemented by theo- 



^ "National Nuclear Data Center, information extracted from the 
NuDat 2 database", http://www.nndc.bnl.gov/nudat2/ 



retical beta-decay rates (IMoUer et al.ll2003h . Neutrino 
absorption on nucleons and nuclei is also taken into ac- 
count. We assume a constant neutrino luminosity for 
neutrinos and anti-neutrinos (Ly^ = ^v, =5x10^^ erg), 
and Fermi-Dirac spectra with a temperature consistent 
with the hydrodynamical simulations. 

In order to understand the nucleosynthesis evolu- 
tion, it is useful to look at the time variations of 
the mean lifetimes for -decays, (n,y), (p,y), (n,p), 
(y, n), and (y, p) reactions. These mean lifetimes (or 
average timescales) are defined as 
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where /1^+(Z,A) denotes the decay rate of nucleus 
(Z,A), {crv}x(Z,A) the reaction rate for process x on 
nucleus (Z, A), and Ayn(Z,A), and Ayp(Z,A) the photo- 
dissociation reaction rates with emission of neutron 
and proton, respectively. In these equations, Na is 
Avogadro number. The average is taken over the heavy 
nuclei (Z > 2): 



Yh= Yj ^(^'^)- 



(7) 



Z>2,A 



Another useful quantity to understand the nucle- 
osynthesis evolution is the reaction flow between two 
nuclei / and j defined as 



F,j^Y(i^j)-Y(j^i). 



(8) 
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The quantity Y(i j) denotes the change in abun- 
dance of nucleus / due to all reactions connecting nu- 
cleus / with nucleus j. The largest reaction flows at 
each time step indicate the key reactions and thus the 
nucleosynthesis path. 

4. Results 

4.1. Constant temperature 

We study here the impact of the wind termination 
varying the temperature at which it occurs. We as- 
sume that the velocity is subsonic and hence there is no 
jump in temperature, density, and velocity. The difl'er- 
ent temperature evolutions and their nucleosynthesis 
are shown in Fig.[2l The original trajectory (solid black 
line) is also included for completeness and comparison 
purposes. In the modified trajectories the wind termi- 
nation is at = 1, 1.5, 2, 2.5, and 3 GK. The resulting 
abundances change significantly with temperature T^t- 
There is an optin ial T^t around 2 GK for producing 
heavier elements dWanajo et al.ll2Ql lab . Moreover, the 
nucleosynthesis evolution is very diff'erent for temper- 
atures higher and lower than this optimal temperature. 
This can be seen in Fig. [3] where the relevant averaged 
timescales (see Sect. [3]) are shown for the trajectories 
with Twt = 1, 2, 3 GK. Initially, there is (p, y) - (y, p) 
equilibrium in all cases. The timescales for the (p, y) 
and (y, p) reactions are much shorter when the wind 
termination occurs at higher temperatures (see bottom 
panel of Figure [3]). 

For the wind termination at low temperatures, 
Twt = 1, 1.5 GK, only elements up to Germanium 
are produced in substantial amounts. In these evo- 
lutions matter expands very fast. Therefore, it only 
stays for a short time in the temperature rage of 
1.5-3 GK where charged-particle reactions can ef- 
fectively synthesize heavier nuclei. In addition, matter 
rapidly reaches large radii where the antineutrino flux 
is rather low. Similar conditions are found in explo- 
sions of low mass progenitors w here the expansion of 
the supernova eject a is very fas t (Hofl^m an et al.ll20Q8l: 
Wanaio et aDl2QQ9l: [Roberts et al. 2QlQ: IWanajo etal ' 



naio 

2011a|) . Notice that the wind termination at a slightly 



higher temperature, 1.5 GK instead of 1 GK, allows to 
reach somewhat higher mass number. 

For the wind termination at high temperatures, one 
would expect that photo-dissociation stops the synthe- 
sis of heavier elements. In this case, the abundance 
would continuously shift towards lower mass number 
as temperature increases. However, in our calculations 
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Fig. 2. — Difl'erent temperature evolutions are shown 
in the upper panel. The solid black line corresponds to 
the original trajectory. The abundances for these evo- 
lutions are presented (with same colors and line styles) 
in the middle and bottom panels versus mass and pro- 
ton numbers, respectively. 



we see an abrupt change in the abundances when the 
wind termination temperature exceeds certain value. 
This points to a key reaction acting as a bottleneck 
at high temperatures. We have identified such a re- 
action using the reaction flows introduced in Sect. [51 
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Fig. 3. — Averaged time scales of the most relevant 
reactions indicated in the labels. The three panels cor- 
respond to the wind termination at 1, 2, and 3 GK from 
the top to the bottom. 



peratures, the reaction ^^Cu(p, Q')^^Ni dominates over 
^^Cu(p, y)^^Zn while at low temperatures the contrary 
is true (see Fig. [5]). When the ^^Cu(p, a) reaction 
dominates, the nucleosynthesis flow is confined into a 
closed NiCu cycle. A similar behavior ha s been found 
in the rp-process for the SnSbTe cycle dSchatz et al 
2QQlh . 
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Fig. 4. — Arrows indicate the flow of difl'erent reac- 
tions normalized to the strongest one. The colors and 
sizes of the arrows are proportional to the flow. Abun- 
dances are also shown with colors and stable nuclei 
are indicated with a dot. The upper panel corresponds 
to the early evolution when the temperature is around 
3.3 GK and in the bottom panel it has dropped down 
to ~ 2.2 GK. 



Figure H] shows these flows at two difl'erent tempera- 
tures for the evolution with T^t = 2 GK. The upper 
panel corresponds to T ^ 3.3 GK, and the bottom 
one to r ^ 2.2 GK, both represent conditions be- 
fore the constant temperature phase. Notice that the 
reactions that determine the nucleosynthesis flow are 
diff'erent at high and low temperatures. At high tem- 



In the calculation with T^t = 3 GK, when tem- 
peratures are still above ~ 3.2 GK (see Fig. [5]) the 
heaviest nucleus produced is ^^Ni due to the NiCu 
cycle. Once the temperature drops, the cycle opens 
and the path reaches ^"^Ge. However, during the high 
constant temperature phase, the triple alpha reaction 
maintains a continuous pr oduction of seed nucl ei from 
fight nuclei with A < 1 dWanajo et aH 12011 ah . This 
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Fig. 5. — Ratio of proton capture rates on ^^Cu, with 
^(P,y) = pYpNA{crv)^p^y) and A^p^a) = pYpNA{crv)(p^a)- 
At high temperatures ^^Cu(p, a)^^Ni dominates while 
for lower temperature ^^Cu(p, y)^^Zn dominates. Dot- 
ted line indicates where rates are equal. 



leads to a reduction of the proton abundance and thus 
of the neutrons produced by antineutrino absorption. 
Notice that the ratio of neutrons produced per seed 
nucleus is a useful guide of the strength of the vp- 
process dPruet et al.l l2006V Both effects (the NiCu cy- 
cle and the reduced neutron-to-seed ratio) result in a 
complete shutdown of y/7-process nucleosynthesis. 

When the wind termination temperature is very 
high (^^ 3 GK) or very low (^^ 1 GK), elements beyond 
Germanium are hardly synthesized. However, the final 
abundances are diff'erent for these two extreme cases. 
When the wind termination takes place at high tem- 
peratures (> 3 GK) matter accumulates in iron group 
nuclei due to the NiCu cycle and to the continuous pro- 
duction of seed nuclei. While for the wind termina- 
tion at low temperatures matter moves fast to large dis- 
tances where the antineutrino flux is not large enough 
to produce a y/?-process. In addition, the low temper- 
atures inhibit the production of seed nuclei. There- 
fore, the final proton abundance is significantly higher 
in the evolution with the wind termination at tempera- 
tures < 2 GK. 



4.2. Reverse shock 

When matter moves supersonically there is a jump 
at the wind termination in temperature and density, 
also called reverse shock. Here, we analyze the im- 
pact of such a jump on the y/7-process nucleosynthe- 
sis using the trajectories shown in Fig. [S] The trajec- 
tory without jump corresponds to the evolution with 



Twt = 2 GK presented in the previous section. The tra- 
jectory with jump is chosen such that the temperature 
increases after the wind termination reaches ^ 2 GK 
(same temperature as in the trajectory without jump). 
The constant temperature phase is the same for both 
evolutions (A^ = 0.5 s). However, the final abundances 
are very diff'erent as shown in Fig.O 
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Fig. 6. — The upper panel shows the temperature 
(thick lines) and density (thin lines) evolutions with 
jump (solid black) and without jump (dashed green) at 
the wind termination. The corresponding abundances 
vs. mass number are presented in the lower panel. 

In the evolution with jump, the expansion contin- 
ues very fast between the moment the temperature 
reaches 2 GK (marked with an A in the temperature 
curve. Fig. O and the position of the wind termina- 
tion shock (marked with B). During this phase, there 
is not enough time for antineutrino absorption on pro- 
tons to produce the necessary amount of neutrons to 
reach heavier nuclei. Already in this phase, between 
A and B, matter starts to beta decay towards stabil- 
ity because temperatures become too low for proton- 
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capture reactions. The temperature rise at the wind 
termination, from point B to C, increases the effective- 
ness of proton-capture reactions. This results in the 
matter flow moving again away from stabiHty and to- 
wards heavier nuclei. Note that the temperature and 
matter distribution after the wind termination are very 
similar for both trajectories (points A and C). 
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Fig. 7. — Same as Fig.[3]for the evolution with a jump 
at the reverse shock. 

The relevant timescales for the trajectory with jump 
are presented in Fig. [T] This can be compared to the 
middle panel in Fig. [3] that corresponds to the trajec- 
tory without jump. After the wind termination, proton- 
captures are faster for the evolution without jump. 
More significant are the difl'erences in (n, y) and (n, p) 
timescales: both are much shorter in the evolution 
without jump (middle panel in Fig. [3]). These difl'er- 
ences in the relevant timescales have a big impact on 
the abundances. 

The reactions involving neutrons, i.e., (n,y) and 
(n,p), depend on the neutron density which is shown 
in Fig. [8] for the two evolutions. In the trajectory with 
jump (solid line), the neutron density remains lower 
than in the trajectory without jump at all times, even at 
the wind termination (feature at ^ 0.02 s). In the re- 
process there is an equilibrium between neutron cap- 
ture and neutron production by antineutrino absorption 
on protons. This implies that 



dYn 
dt 



Here Ay^ is the electron antineutrino absorption rate 
and <(Tv)(z,A) is the sum of reaction rates for (n, y) and 
(n, p) reactions for nucleus (Z, A). Therefore, the neu- 



tron density in equilibrium is given by 



(10) 



The nucleosynthesis path is very similar for both tra- 
jectories considered here (evolution with and without 
jump). Therefore, only small variations are expected 
in the denominator. Consequently, the difl'erence in 
the neutron densities (see Fig. [8]) is due to the neutron 
production by antineutrino s. The deceleration of the 
expansion at the wind termination occurs at a smaller 
radii for the trajectory without jump (see bottom panel 
of Fig. [8]). The production of neutrons is hence less 
eflScient in the case with jump because matter reaches 
larger radii where the neutrino flux is reduced due to its 
dependency. The difl'erence in the abundance pat- 
tern for nuclei above mass number A 80 (see Fig.[6l) 
is due to this difl'erence in the efficiency of neutron pro- 
duction. 

4.3. Late-time evolution 

We explore the impact of the dynamical evolution 
after the wind termination. In this section, we assume 
that the wind termination occurs at T^t = 2.0 GK and 
we vary the parameter Ar. The quantity Ar character- 
izes the timescale for the transition from an expansion 
with constant temperature and density (during which 
V oc r"^) to a constant velocity expansion with p oc r~^. 
We choose the values At = 0.0, 0.5, and 1.0 s, moti- 
vated by the anisotropic evol ution of the ejecta in 2p 
hydrodynamical simulations ( Arcones & Jankall201 1 ). 
There is a smooth transition between the two extreme 
cases of A^ = 1.0 s and At = 0.0 s, which can be 
seen in the intermediate case of A^ = 0.5 s. The latter 
(At = 0.5 s) was used in previous sections. 

The diff'erent temperature evolutions (top panel) 
and the resulting nucleosynthesis (middle and bottom 
panels) are shown in Fig.[9l While all three cases pro- 
duce similar abundance distributions, the details de- 
pend critically on the value of At. 

For a relatively long phase of constant temperature 
(At = 1.0 s; left column of Fig.[TO|), the (p,y) - (y,p) 
equilibrium lasts for quite a long time and allows for 
many (n, p) and (p, y) reactions to occur. This drives 
the matter to higher mass number (up to A ^ 110), 
where matter accumulates in Sn (Z = 50 closed shell). 
Due to the extended period of (/?, y) - (y, p) equilib- 
rium, significant abundances of nuclei with 64 < A < 
110 are synthesized (see line C in Fig. [TOl left col- 



7 





10 15 20 25 30 35 40 45 50 55 60 



Fig. 10. — The left column is for the case of ISX = 1.0 s, the right column is for = 0.0 s. The top row shows average 
timescales for different reaction types as function of time. Abundance distributions at selected times are shown as 
function of mass number (middle row) and as function of atomic number (bottom row) for both cases of A^. For the 
case A^ = 1.0 s, the abundances are shown at T = 3.0 (A), 2.0 (B,C), 1.6 (D),1.0 GK (E) and also the final abundances 
(F). The abundances at the beginning (B) and at the end (C) of the constant temperature phase diff'er quite significantly. 
For the case ^t = 0.0 s, the abundances are shown at T = 3.0 (A), 2.0 (B), 1.9 (C), 1.5 (D), 1.0 GK (E) and also the 
final abundances (F). Note that while the abundances at (C) for A^ = 1.0 s and at (C) for A^ = 0.0 s are similar, the 
temperature is slightly lower for (C) because there is no constant temperature phase in this case. 
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time [s] 



Fig. 8. — Evolution of neutron density (upper panel) 
and radius (bottom panel) for the trajectories presented 
in Fig.[6l The solid black line corresponds to the evo- 
lution with jump, the dashed green line to the evolution 
without jump. 

umn). Also during this phase there is a continues pro- 
duction of seed nuclei at the expenses of alpha parti- 
cle (see Fig. 11). Once the temperature drops below 
1.5 GK (line D), the production of seed stops. In 
addition, proton-capture reactions become hindered by 
the Coulomb barrier and -decays become faster than 
(n,p) reactions (see upper left panel in Fig [9]). The 
nucleosynthesis ceases to efficiently proceed to higher 
masses. During this phase, matter starts to decay to- 
wards stability. A small number of late time neutron 
capture reactions smooths the abundance distribution 
in the mass range 80 < A < 1 10 (line E). Once the tem- 
perature drops below 1.0 GK (lines E and F), the abun- 
dances as function of mass number do not change any- 
more, as at this time -decays dominate. The closed 
shell at Z = 50 acts as barrier for the nucleosynthe- 
sis. This can be seen in the enhanced abundances at 




z 



Fig. 9. — Different temperature evolutions after the 
wind termination shock which is at 2.0 GK are shown 
in the top panel. Values for At are 0.0 s (solid black), 
0.5 s (dashed green), and 1.0 s (dotted red). Abun- 
dances for these evolutions are shown versus mass 
number A (middle panel) and versus atomic number 
Z (bottom panel). 

A = 110 (originating from ii0'111'1123j^^^ ^^^i 
abundance distribution. 

For the case = 0.0 s (right column in Fig. [TOb . 
the situation is very different. The initial phase of 
(p, y) - (y, p) equilibrium is much shorter than in the 
case of At = 1.0 s. In this case nuclei up to mass 
number A ^ 110 are synthesized (line C') the tem- 
perature has already dropped below 2 GK. Due to the 
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faster temperature decline, the production of seed nu- 
clei is less efficient than in the case of = 1.0 s (see 
Fig.[TT]). This leads to lower abundances of nuclei such 
C, 1^0, 20Ne, 28si, and ^OCa for = s (see 



as 



12, 



Fig.O. The lower abundance of seed nuclei leads to a 
higher neutron density (see Eq.[TO|) as in both cases the 
neutron production rate ( Ay^Yp) is similar. Therefore, 
in the evolution with A^ = s, there are more neutrons 
available per seed nuclei (Fig. [TT]). The neutron cap- 
ture reactions, i.e. (n,y) and (n,p), become faster than 
(p, y) reactions at temperatures of T < 1.5 GK. These 
neutron captures move matter to the neutron-rich side 
of stability. This, together with the earlier less efficient 
production of seeds, results in a depletion of the abun- 
dances in the region of 64 < A < 1 10 (compare black 
solid and red dashed line in Fig. [9]). Towards the end of 
the evolution, (n, y) reactions are the fastest ones, even 
faster than -decays. During this phase, the atomic 
number Z remains unchanged while the mass number 
increases (lines D, E, and F in Fig.[TO|). The signature 
of matter moving to stability via (n, y) reactions can 
be seen in the increased abundances at A = 124, cor- 
responding to ^^^Sn, the most neutron-rich stable Sn 
isotope. 

5. Conclusions 

We have studied the impact of the supernova dy- 
namical evolution on the y/7-process and how the wind 
termination affects the synthesis of elements beyond 
iron. For a robust and strong production of nuclei with 
64 < A < 110 there is an optimal win d termination 



tempera ture at ^^2 GK as it was found by lWanajo et al 



dZOlla"). Also in agreement with their work, we have 
shown that if the wind termination occurs at low tem- 
peratures (r^t ~ 1 GK), matter stays too short time 
in the optimal temperature range for y/7-process nucle- 
osynthesis, ~ 1-2 GK. Moreover, the electron antineu- 
trino flux quickly becomes too small to produce the 
necessary neutrons to overcome the -decay waiting 
points. This hinders the y/?-process and consequently 
the efficient synthesis of heavy nuclei. Therefore, the 
wind termination temperature determines the heaviest 
elements produced. 

We have identified an end point cycle that is key 
at high temperatures (around 3 GK). The close NiCu 
cycle inhibits the production of elements heavier than 
56Ni. The reactions in this closed NiCu cycle are the 
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Fig. 11. — Abundances of a particles (top panel), 
abundances of heavy nuclei with Z > 2 (middle panel), 
and neutron densities (bottom panel) for Ar = 0.0 s 
(solid black), 0.5 s (dashed green), and 1.0 s (dotted 
red). 



following: 



56^1 



Ni(n, prCoip, yy'Niin, pf'Co(p, yf^'Ni 
^^Ni(p,yf^Cu(p,af^Ni. 

At high temperatures the reaction ^^Cu(p, a) ^^Ni pre- 
vents the synthesis of elements beyond the iron group. 
When the temperature drops below ^ 3 GK, ^^Cu(p, y) 
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Zn becomes more effective than Cu(p,a) Ni. 
This leads to breakout from the NiCu cycle and the 
flow of matter can continue towards heavier nuclei. 
The temperature dependence of these two reactions is 
critical because it sets the temperature at which the 
y/7-process can start to synthesize elements beyond 
Nickel. If this cross-over temperature were slightly 
higher, matter would be closer to the proto-neutron 
star and thus under higher neutrino flux when the path 
starts to move towards heavier nuclei. This will sig- 
nificantly increase the efficiency of the y/7-process. 
Therefore, our results clearly motivate further investi- 
gation of these two key reactions: ^^Cu(p, a) ^^Ni and 
^^Cu(p, y)^^Zn. Particularly relevant are the branching 
ratios for the decay by alpha and gamma emission of 
compound states in ^^Zn above the proton separation 
energy. 

We have also explored for the first time the im- 
pact of the dynamical evolution after the wind termi- 
nation on the synthesis of heavy elements by the vp- 
process. Hydrodynamical simulations show that after 
the wind termination there is a transition from an ex- 
pansion with almost constant temperature and density 
to a phase with almost constant velocity. The duration 
of the constant temperature phase. At, strongly aff'ects 
the final abundances of heavy nuclei. We have found 
that for = 0.0 s the final flow of matter to stabil- 
ity occurs by neutron capture reactions, i.e (n,y) and 
(n, p), and not by -decays. When the phase of con- 
stant temperature is very short (A^ < 0.5 s), the abun- 
dances of nuclei with 64 < A < 110 are lower. This 
leads to significantly higher neutron densities at later 
times. Depending on the neutron density the matter 
moves to stability either by -decays (At = 1.0 s) 
or by neutron captures. These two different ways of 
reaching stability leave a distinct fingerprint in the fi- 
nal abundances. 

We have investigated in detail the impact of dynam- 
ical evolution on the y/7-process nucleosynthesis using 
individual characteristic trajectories. In order to pre- 
dict the complete y/7-process yields from a supernova 
simulation, one will need to integrate over all proton- 
rich ejecta. 

In summary, the supernova dynamics as well as in- 
dividual reactions such as the proton capture reactions 
on ^^Cu determine how high in proton number the vp- 
process can proceed. The dynamical evolution just af- 
ter the wind termination can greatly affect the nucle- 
osynthesis evolution towards stability as it determines 
the late neutron density. Our results provide a link 



between nucleosynthesis in proton-rich winds and the 
dynamical evolution of the ejected matter and motivate 
further theoretical and experimental eflfort on under- 
standing key reactions. 
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